The Compressible Laminar Boundary Layer 
on a Yawed Infinite Wing* 


L. F. CRABTREE, B.Sc., Grad.R.Ae.S. 


Summary: A study, of the laminar boundary layer on a yawed infinite wing in a 
compressible fluid is made. A method is presented for the investigation of compressi- 
bility effects by direct solution of the linearised equations of motion, first considerably 
simplified by the extension of a transformation due to Illingworth and Stewartson, 
which assumes a heat-insulated surface and a fluid of unit Prandtl number and a 
coefficient of viscosity which is proportional to the absolute temperature. As an example 
of this procedure the boundary layer near a stagnation point is calculated. The 
simplifications necessary for an extension of the momentum method, which is available 
for a yawed infinite wing in an incompressible fluid, are discussed, and this method is 

recommended for a first evaluation of the effect of compressibility. 


Notation 
a_ velocity of sound 
c a constant in equation (8) 
d a constant in equation (18) 
f Blasius chordwise-flow function (equation (9) ) 
g _ velocity profile of the spanwise boundary layer (equation (9) ) 
K a constant in equation (21) 


I non-dimensional shearing-stress parameter (see equations (30) 
and (34) ) 


_ _ aw (6) 
_M_ Mach number 


T absolute temperature 


= modified Pohlhausen parameter (equation (29) ) 


u_ chordwise velocity component inside the boundary layer 

U ___chordwise velocity component outside the boundary layer 

v spanwise velocity component inside the boundary layer 

V spanwise velocity component outside the boundary layer 
co-ordinate measured along the surface in the chordwise direction 


x 
y co-ordinate measured in the spanwise direction 
Zz 


co-ordinate measured normal to the surface 


*Based on part of a thesis presented in partial fulfilment of the requirements for the degree 
of Ph.D., Cornell University. 
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ratio of specific heats 
a Blasius variable (equation (9) ) 


(4 4) dz=momentum thickness of the boundary layer in 
p,U U the chordwise direction 


sea (1 - 7 dz=momentum thickness of the boundary layer in 
pV the spanwise direction 


5 (1- A dz=“mixed” momentum thickness of the boundary 
layer 


coefficient of viscosity 
kinematic viscosity 
density 
4(y-1) My’ 
= x/c (equation (28) ) 
shearing stress at the wall 
stream function (equations (6) and (9) ) 


Subscripts 


refers to stagnation conditions 

refers to conditions at the edge of the boundary layer 
refers to a related incompressible flow 

refers to a chordwise component 


refers to a spanwise component 


1. Equations of Motion 


For incompressible flow over a yawed infinite wing, the chordwise flow, both 
inside and outside the boundary layer, is independent of the spanwise flow, as has 
been shown by Prandtl’, Jones’, and Sears’. For the same wing in compressible 
flow, however, Struminsky“’ and Moore’? have shown that this “ independence 
principle” does not hold (with the exception of isothermal fluids). This may be 
seen physically, since the temperature affects the chordwise velocity profile by 
changing the density, and the temperature itself is dependent upon the resultant 
velocity, i.e. upon both the spanwise and chordwise velocity components. 


By the use of a transformation of variables and of the velocity distribution 
outside the boundary layer, due to Illingworth" and Stewartson"?’, the equations 
of motion for the two-dimensional compressible boundary layer problem may be 
reduced exactly to those of a related incompressible flow (see also Ref. 1). To 
obtain this result it is necessary to prescribe a heat-insulated surface and a gas 
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of Prandtl number unity whose coefficient of viscosity varies directly with the 
absolute temperature. It is evident, by the reasoning already applied to the question 
of existence of an independence principle, that the Illingworth-Stewartson trans- 
formation will not reduce the boundary layer equations of motion for compressible 
flow over a yawed infinite wing exactly to an incompressible case. However, since 
great simplification of the two-dimensional problem was achieved by this trans- 
formation, it will be used again here, and at least some advantage may be expected. 
It will be necessary to prescribe the same restrictions of a heat-insulated wall and 
a gas of Prandtl number unity and whose viscosity is proportional to the absolute 
temperature. 


The properly extended transformation of Illingworth and Stewartson is applied 
to the boundary layer equations of motion of a yawed infinite wing in the Appendix. 
New co-ordinates x; and z; are introduced according to the equations 


3y-1 


(1) 


(2) 


where x is the co-ordinate measured along the surface in the chordwise direction 
and z is measured normal to the surface. The density p and speed of sound a are 
characterised by subscripts 0 referring to stagnation conditions, or 1 referring to 
conditions at the edge of the boundary layer. Absence of a subscript indicates 
ambient conditions inside the boundary layer. 


The chordwise component U of the velocity at the edge of the boundary layer 
is transformed by writing 


Equations (1) to (3) reduce the chordwise and spanwise equations of motion to 
"0 Oz, Ox0Z, Qx,0z° 2 V ] @) 


respectively, where the stream function y¥ is defined by 


pu=p, 
pw=—p,dp/ax 


My*=V? | ve) 
87 


|_| 
0 


L. F. CRABTREE 


i.e. My is the Mach number of the spanwise flow at the edge of the boundary layer 
with respect to the speed of sound at a stagnation point of the chordwise flow. 


Except for the last term on the right hand side of equation (4) it is seen that 
this equation and equation (5) are the equations of motion for a related incom- 
pressible flow. Further, the parameter 


-1 

is independent of the chordwise co-ordinate x, since My is constant, by virtue of 
equation (7), if V, the spanwise velocity at the edge of the boundary layer, 
is constant. 


We consider it as an extreme case when the Mach number of the spanwise flow 
component reaches unity, since our treatment cannot take into account shock wave-— 
boundary layer interaction. With this restriction the maximum possible value of 
w is 0-2, for air with y=1:4, and M,=1. 


This suggests a linearisation of the problem by a series expansion in powers 
of the parameter w. 


It must be emphasised at this point that the additional term involving © in 
equation (4), the chordwise equation of motion, does not represent the complete 
effect of the spanwise velocity on the chordwise velocity profile; a part of the 
compressibility effect of the spanwise on the chordwise velocity is implicit in the 
remaining terms of the equation. This state of affairs stems from the use of the 
extended Illingworth-Stewartson transformation, as the ratios p/p, and T/T, which 
are involved are calculated from the resultant velocities. 


2. Flow Near a Stagnation Point 


The situation was investigated for the case where the related incompressible 
flow is that of flow near a stagnation point, i.e. 


By utilising a transformation of variables due to Blasius (Ref. 8, p. 139) 


v= x f () 


v=Veg(n) 


where V is, of course, constant, equations (4) and (5) reduce to 


{1+0[1-(3) ]} . (10) 
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with boundary conditions f(0)=f(0)=0, f(co)=1, 
and ge” +fe’=0 ‘ 
with boundary conditions ¢(0)=0, g(oo)=1. 
In these equations primes denote differentiation with respect to ». 
We now expand the functions f and g in series of powers of w:— 


f= f, + af,+..... (12) 


f, is found to satisfy the equation 


as given by Blasius and Hiemenz. The solution may be found tabulated in Ref. 9. 


Thus f, is the same function as in the incompressible case, and if first order 
terms in © are neglected, the chordwise flow function of our transformed system 
is independent of the spanwise flow function g. 


The equation for f, is 
Also, g, and g, satisfy the equations 


respectively. 


Equation (16), together with its boundary conditions g,(0)=0 and g, (00)=1, 


is exactly that given and solved by Sears’ and Cooke”. Its solution is 
| 
0 0 


with the constant d=0-570, and has been tabulated by Cooke”. 


With the known values of f, and g,, equation (15) may be integrated 
numerically for f,. This has been done by Adams’s method and the results are 
given in Table I. The effect of the first-order correction term wf, on the two- 
dimensional solution f, is shown in Fig. 1 for »=0-2. The corresponding effect 
on the first derivative of f, (or the effect on the velocity) is plotted in Fig. 2. It is 
seen that the velocity correction term wf,’ is small, 
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Fig. 1. 


Chordwise flow function for stagnation point flow, 
U,=cx,, with w=0°2. 


These results have been transformed back into the physical plane to give the 
actual velocity profiles shown in Fig. 3. The use of the Illingworth-Stewartson 
transformation even at the zero approximation in w (f=f,) is seen to give a signifi- 
cant improvement over the incompressible assumption for the chordwise velocity 
profile. The first approximation, f=f,+f, changes the velocity profile further in 
the same sense. 


Now the Solution of equation (11) is (18) if f, is replaced by f,+f,, with the 
constant d to be evaluated from the boundary condition g(0%o)=1. This numerical 
quadrature has been carried out (see Table I) and the constant is found to be 


for d=0-2. 


Thus the function g, may be calculated without directly solving equation (17), 
and the solution may be written 


" " " 


wg, =0°525 exp ( - | (f, + of,) dn ) dy —0°570 {exp ( fd (20) 


0 0 0 0 
for ©=0-2, neglecting terms of order w’. 


The second term on the right hand side of equation (20) is that tabulated by 
Cooke"*’ as G(Y). An examination of Table I reveals that the effect of the 
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Fig. 
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U,=cx,, with @=0°2. 


Chordwise velocity function for stagnation point flow, 


CASE 


APPROXIMATION~_| 
ZERO ORDER APPROXIMATION 


INCOMPRESSIBLE 


Chordwise velocity profiles for stagnation point flow, U;=cx,, according 
to the incompressible assumption, and also the zero and first order 
- approximations to the compressible case, with « =0-2. 
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TABLE I 
CHORDWISE AND SPANWISE PROFILES (STAGNATION POINT FLOW) WITH o=0-2 


f,+of, 8, g, 


0 0 0 
0-0060 0:0041 0:0811 0:0525 
0:0233 0:0162 0-1582 0-1052 
0-0510 0:0357 0:2313 0-1575 
0-0881 0:0625 0:3005 0-2098 
0:1336 0:0958 03657 0-2617 
0° 1867 0°1355 0°4271 0°3131 
0:2466 01812 0:4847 0:3638 
0°3124 0:2324 0:5386 04136 
0-3835 0-2888 0:5887 0°4622 
0°4592 0:3501 0°6352 05095 
0°6220 0°4856 0:7176 0:5988 
0-°7966 0:°6362 0°7864 0:6798 
0:9798 0-:7994 0°8425 0:7511 
1:1688 0-9724 0-8870 0-8120 
1:3619 1-1533 09211 08621 
1:5577 1:3402 0:9469 0:9021 
1-7552 1:5316 0:9651 0:9327 
1-9537 1:7259 0:9777 0-9552 
2°1529 1:9224 0:9865 0:9713 
2°3525 2°1203 09921 0-9822 
2°5522 2:3190 0:9955 09894 
2°7521 2°5183 0-9975 0:9939 
2-9520 2:7179 0:9987 0:9967 
3°1518 2°9175 0-9993 0-9983 
3-3518 1-0000 3-1175 0:9997 0:9992 
3°5518 1-0000 3-3174 0-9999 0-9996 
3°7518 1-0000 3°5174 0-99995 0-9999 
3-9518 1-0000 3-7174 1-0000 1-0000 1-0000 
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first-order velocity correction term, wg,, is small. Thus the “universal” spanwise 
velocity profile") is also found for the compressible case, to a good approximation. 


Having obtained a first-order approximation to the solution for f, the equations 
(10) and (11) may now be solved directly by the standard trial and error method. 
Several numerical integrations to infinity will be required and the value of f” (0) 
is adjusted after each integration until the boundary condition at infinity is satisfied. 


3. The General Falkner-Skan Case 


The case of flow near a stagnation point is a special case with K=1 of the 
more general solution given by Falkner and Skan for the two-dimensional 
incompressible flow 


(see Ref. 8, p. 140). 
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The equations for steady flow in the boundary layer, (4) and (5), are satisfied 
by putting 


+ U; 


With these substitutions equation (4) becomes a differential equation for f, 


{1+0[1- (3) |} 
which may be compared with the equation obtained by Falkner and Skan 


The equation of motion for the spanwise flow, (5), may be reduced by the 
same substitutions together with 


v=Vg(n)_ . . (26) 
and the result is +4(K+1)fe=0,. ‘ ‘ ‘ 


which is exactly the result obtained by Cooke" in the incompressible case. 


The linearisation procedure adopted for the special case of stagnation point 
flow in Section 2 may, therefore, be utilised again here, and it seems evident that 
the compressibility corrections in this general case will have the same order of 
magnitude as in the special case. 


In principle then, the solution of the laminar boundary layer on a yawed wing 
in a compressible fluid may be obtained by a direct solution of the differential 
equations of motion, as has been indicated. But generally this method is too 
complicated, so that it seems desirable to use approximate “momentum” or 
“Pohlhausen” methods for the boundary layer development. 


4. The Use of Momentum Methods 


4.1. The isothermal compressible case 


If the fluid is assumed to be isothermal the independence principle holds, as 
has been noted previously. Thus the momentum methods developed in Ref. 1 for 
the incompressible case may be directly extended to the case of a compressible fluid. 


The isothermal assumption is, however, a further approximation, and its effect 
was investigated by calculating the chordwise characteristics of the boundary layer 
at whose edge the velocity distribution is given by 


U=U, (€-&) iw’ % (28) 
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INCOMPRESSIBLE 
COMPRESSIBLE 
COMPRESSIBLE 


Momentum thickness of the laminar boundary layer for the velocity 
distribution U=U,, 


representing the flow over a certain class of cylinders, where £=x/c, c being a 
typical dimension. The results presented in Ref. 1 are used to calculate the com- 
pressible case, and the same results are easily modified to include the additional 
assumption of an isothermal fluid. The relevant formulae for the limiting case, 
y=1, are 


2 
6,7=0-45 v, exp(M’) U-* exp ( U*dx . (29) 


(30) 


From the values of m, given by equation (29) the corresponding values of /; 
are read from the universal curve of Ref. 1. The shear stress at the wall may then 
be calculated from 


The values of the momentum thickness @,, and of -, the shear stress at the 
wall, for the compressible and isothermal compressible cases are plotted in Figs. 4 
and 5, in comparison with the incompressible results. 


This comparison shows that the isothermal compressible results are closer to 
the actual compressible values than are the incompressible results. 
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Wall shear stress for the velocity distribution U =U, (€—&°). 


4.2. Non-isothermal fluid—further simplifications 


The momentum method may be extended, however, to the case of a non- 
isothermal fluid, so that although the complete dependence of the chordwise upon 
the spanwise flow, in the compressible case, is not represented by the method, at 
least a part of the dependence is accounted for in the solution. This part dependence 
is obtained by making use of the Illingworth-Stewartson transformation, as 
explained in Section 1 and, provided that the first order term in involving the 
spanwise velocity profile is neglected, the equation of motion for the chordwise 
flow is reduced to that of a related incompressible flow. 


The results obtained in Part 4 of Ref. 1 for the two-dimensional compressible 
case may therefore be taken over completely for the calculation of the chordwise 
characteristics of the boundary layer in the present problem. 


The calculation of the spanwise flow in the compressible case requires some 
further consideration. It may be shown that by use of the extended Illingworth- 
Stewartson transformation, the spanwise momentum equation reduces to 


This represents the spanwise momentum equation for a possible incompressible flow 
pattern over a yawed infinite cylinder, and it is therefore reasonable to invoke the 
aid of the assumptions made for such a case in Ref. 1. Thus, having calculated 
6, and m, 9, (the momentum thickness in the spanwise direction for the related 
incompressible flow) may be calculated by any of the methods presented therein. 
The universal spanwise profile then gives v/V directly, and the remaining properties 
of the boundary-layer may be calculated from the following formulae :— 
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y+} 
2ly—1) 
0,=6, (2°) 
1 
dv dv 
1 (dv l(m=0) 
or, since (35) 
and l(m=0)=0-225, 
Sy—1 
2y—1 
vi 0 
or r4=0°225 (72). 
v 0 


Thus, the compressibility effect for yawed wings can partially be taken into account 
by the momentum method without any essential complication. 


Conclusions 


In principle, the solution of the laminar boundary layer on a yawed infinite 
wing in a compressible fluid may be obtained by a direct solution of the linearised 
differential equations of motion. These equations are first considerably simplified 
by the extension of a transformation due to Illingworth and Stewartson which 
entails the assumptions of a heat-insulated surface and a fluid of unit Prandtl 
number and a coefficient of viscosity which is proportional to the absolute 
temperature. Even so, this method is generally too complicated, and therefore the 
use of an approximate momentum method is considered for the boundary layer 
development. 


First, the further assumption of an isothermal fluid allows of the direct 
extension of the momentum method which is available for incompressible flow. 


Secondly, the linearisation procedure is used to extend this method to the case 
of a non-isothermal fluid. The momentum method may thus be used to account 
partially for the compressibility effect on yawed wings without any essential 
complication. 
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APPENDIX 


TRANSFORMATION OF THE BOUNDARY LAYER EQUATIONS OF A 
YAWED INFINITE WING IN A COMPRESSIBLE FLUID 


Using the co-ordinate system defined in Section 1, the equations of motion for 
steady flow are 


together with the equation of continuity 


The equation expressing conservation of energy is 


for a Prandtl number of unity. 
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If there is no heat transfer between the gas and the wing surface, the simple 
integral of equation (42) may be used, 


2 


The Illingworth-Stewartson transformation is‘ 


Note that 
where V’=the resultant velocity at the edge of the boundary layer. 


Equation (41) implies the existence of a stream function ¥ such that 


Assuming further that . (46) 


i.e. the coefficient of viscosity is directly proportional to the absolute temperature, 
and substituting from the preceding equations into equation (38), gives 


2 ! 2 92 3 93 
a, & | + y-1 dx __ oy (47) 


2y 


so that equation (47) becomes 


ox \a,/ p, p, dx 1+ 
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Therefore 


p, dx 


Oz, | Ox; Oz? + (U?+V?—v?) | (50) 


2a, 


Now U\=(a,/a,)U, from equations (44). 


Therefore dU; _ 


and _y, ay dp, @ 


dx, a, p, dx 


since for V’? =U? + V? with constant V, then V’ dV’ /dx,=U dU /dx, and substitution 
into the right hand side of equation (50) gives 


Oy dy ow | y-1 


where M,?=V?/(a,?—4(y-1) V’), 


i.e. My is the Mach number of the spanwise flow at the edge of the boundary layer 
with respect to the speed of sound at a stagnation point of the chordwise flow. 


Similarly, by substitution into equation (39), 


a, Ox 0% dz p 0z » 0% 
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a, av dv a, dv (2 
Therefore a, 0x a Ox 
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hence, cancellation of like terms in equation (53) gives 
sy-1 
av av y-1 
and by transformation of the co-ordinate x by equations (44), 


Oyov dv 

The equations of motion corresponding to (52) and (55) in the incompressible 
case are 


Ox dz 
ow dv 


(56) 


where u=0y/0z and w= — 0x. 


| = 
Oz Ox Ox Oz? 
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P. J. PALMER, Ph.D., M.Sc., D.LC., A.F.R.Ae.S. 
(Department of Civil Engineering, University of Birmingham) 


Summary: The natural frequency of vibration of curved rectangular plates, with 
both simply-supported and fixed edges, is evaluated for the fundamental extensional 
mode. This mode is thought to be applicable when the plates are excited by a uniformly 
distributed pressure as may occur with a dynamic pressure wave. The natural frequen- 
cies corresponding to this mode increase fairly rapidly with the curvature of the plate. 


The natural frequency of vibration of the plates when the mode is the fundamental 
inextensional mode is also considered. The frequency of this mode is higher than that 
of the extensional mode for small curvatures, but the inextensional mode frequency falls 
slowly with increase in curvature of the plate. Thus, if the curvature of the plate is 
sufficiently large, the frequency of the fundamental inextensional mode is lower than 
that of the extensional mode. 


All the frequencies deduced are found by means of the Rayleigh method, in so 
far as they correspond to initially assumed modes. 


Introduction 

In a recent investigation into the effect of dynamic pressure waves on stressed 
skin structures, it was found necessary to evaluate the natural frequency of vibration 
of curved rectangular plates. Very little information appears to have been published 
on this subject, although as early as 1888 Lamb"? considered the vibration of curved 
bars. In 1928, Den Hartog'?’ evaluated the natural frequency of vibration of 
circular arcs when vibrating in both extensional and inextensional modes. The 
present paper sets out to deal with the natural frequency of vibration of curved 
rectangular plates while vibrating in fundamental modes. 


The natural frequency of vibration of any curved rectangular plate will, of 
course, depend upon the mode into which it deforms. When designing against 
dynamic pressure waves, it is thought that, because the load distribution which 
excites the mode will be approximately uniformly distributed over the plate, the 
mode produced will be the fundamental symmetrical one in which extension of the 
middle plane of the plate occurs. Hence, emphasis is placed here upon this type 
of extensional mode, although an inextensional mode is also investigated. It is 
quite possible that after the initial excitation, by the pressure wave, the plate may 
change its mode of vibration, if one of lower frequency exists, but as regards the 
structural design of the plate and neighbouring structure, it is the mode occurring 
during the first cycle which is relevant. 


Throughout this paper the natural frequencies deduced are those obtained by 
the Rayleigh solution in so far as they correspond to initially assumed modes, so 
that the true natural frequencies may be slightly lower than those deduced here. 


Originally received April 1953. 
[The Aeronautical Quarterly, Vol. V, July 1954] 
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Fig. 1. 
Geometry for curved plate. 
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 Poisson’s ratio 
Strain components in plate 


stress components in plate 


2. The Extensional Mode 


This section is concerned with the natural frequency of vibration of a curved 
plate when the mode is the fundamental symmetrical one, in which extension of 
the centre plane of the plate takes place. 


The case of a flat rectangular plate with simply-supported edges is a classical 
case and appears in textbooks’. The case of a flat plate with fixed edges (fixed 
as regards position and slope) is not so simple, but the solution using the Rayleigh 
method has been worked out by Redshaw™’. The natural frequency of curved 
rectangular plates with edges either simply-supported or fixed can be evaluated, 
using the Rayleigh method, as follows. 


For a curved rectangular plate, as shown in Fig. 1, with axes in the plate of 
Aé and Bo, where A and B are the radii of curvature, the expressions for the strains 
are, from Timoshenko”, 


Add’ A 
“* Boo B 
= 
Bag Aad 


e 


where u, v and w are the displacements in the A#, Bo, and z directions. 


The change of curvature in the A@ direction is 


(2u_ 
A2\00 00? 


(2) 


For the type of mode being considered, extension of the centre plane of the 
plate takes place, and points on the centre plane move radially such that u and v 
are zero. Hence the first term in the expression for the direct strains can be 
expressed as the distance from the centre of the plate multiplied by the change in 
curvature. Also, provided that the thickness of the plate is small compared with 
the radius of curvature, the first term in the expression for change of curvature can 
be neglected compared with the second. 


Thus the expressions for curvature changes are dependent upon w alone, and 
hence the strains can be expressed in terms of w. This leads to the following 
expressions for the strains; — 
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A? 002 

0*w 
B? Og? 

__ 2% ow 

~ AB 0000 


where z is the axis normal to the surface of the plate measured from the 
middle plane: 


In the Rayleigh method, an assumption is made concerning the shape of 
normal deflection of the plate (i.e. the mode). The form taken here for the case 
of simply-supported edges is 


w=Jsin sin : i (4) 


where a and b are the lengths of the sides of the plate. 


This mode is the simplest form which is symmetrical and which satisfies the 
boundary conditions 


=0 at and a/A 
(5) 


“=0 at and 5/B 


It is also that mode which does in fact lead to the exact solution in the case 
of flat plates. 


With fixed or built-in edge supports, the assumed mode will be 
w= J 0°97 (a— (b— Bo)? (6) 


which is a symmetrical mode, which reduces to that used by Redshaw™? in the case 
of flat plates, and which satisfies the boundary conditions given by 


Ow 


at #@=0 and a/A 
(7) 


ow 
w= at and b/B 


The stresses can be expressed in terms of the strains by the following 
relationships. 


+ 


(1-4) 


E 
+ 


G69. 
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12 


afxio> 


n.|sec. 10 


hla =0:025 
h/e-= 9-920. 


h/a = 


0-10 a/y 
Fig. 2. 


Natural frequency of curved square plates. Equal curvature in both directions. 
Simply-supported edges. 


0-10 


Fig. 3. 


Natural frequency of curved square plates. Equal curvature in both directions. 
Fixed edges, 
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| Na = 0-025 
= 0-020 
| Mg = 9-010 


fo] 0-05 O15 
Fig. 4. 


Natural frequency of curved square plates. Curved in one direction only. 
Simply-supported edges. 


The potential energy stored in the plate can be obtained from the expression 


and the kinetic energy, at any time, is given by 
dT =4phw* AB dé do (10) 


The sum of the potential energy and the kinetic energy will be constant, so 
that if the total energy is differentiated with respect to time, a differential equation 
of the following form results. 


The term J is the amplitude of vibration and, as such, will vary with time, whereas 
C is independent of time. From this differential equation an expression for the 


natural frequency of vibration can be deduced. For simply-supported edges, this 
leads to 


and for fixed pee 


These expressions agree with the existing ones for flat plates and it should be 


noted that, in the limit for large curvatures, the type of edge fixing becomes 
immaterial. 
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= 0-025 


Nia = 0-020 


| 


0-008 


0-05 0-10 0-15 


Fig. 5. 
Natural frequency of curved square plates. Curved in one direction only. 
Fixed edges. 


The increase of the natural frequency of a plate with curvature is small for 
very small curvatures, but the increase soon becomes quite large as the curvature 
increases. Typical results are given for light alloy plates in Figs. 2 to 7. These 
results are for square plates with equal curvature in both directions, square plates 
with curvature in one direction alone, and infinitely long rectangular plates curved 
across the width, all three cases being given for both types of edge support. Here 
E has been taken as 10 x 10° Ib./in.*, v as 0-3, and p as 0-1/g slugs/in.* 


3. The Inextensional Mode 


As already mentioned, the natural frequency of the fundamental symmetrical 
mode, in which extension takes place, is simply that mode which is thought to be 
the most relevant when designing against dynamic pressure waves. Other modes 
may occur, notably inextensional modes, in which the natural frequency may be 
less than that of the extensional.mode considered. Here, the fundamental mode 
cannot be symmetrical, and the fundamental asymmetrical mode, in which no 
extension of the middle plane of the plate occurs, will be considered. Whereas in 
the extensional mode the lateral displacements u and v are zero, since the points on 
the centre plane of the plate deform along radial lines, in the inextensional mode 
these displacements will now not be zero, but will be such as to cause no extension 
of the middle plane of the plate. 


Thus, the three strain components will be zero in the middle plane of the 
plate and for the direct strain in the A@ direction this leads to 


ou 
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0-05 0-10 0-15 ol, 
Fig. 6. 
Natural frequency of infinitely long rectangular plates, curved across the width. 
Simply-supported edges. 


The change of curvature in the centre plane of the plate, in the A@ direction, 
then becomes 


1 ) 


and the expression for the corresponding strain is 


z 
62 +w 


0? 
and similarly - +w) 


10 


8 


hla = 0-025 
0.020 
0-018 
Na = 0-010 


h = 


0°05 0-10 0-15 als 
Fig. 7. 


Natural frequency of infinitely long rectangular plates, curved across the width, 
Fixed edges. 
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EXTENSIONAL MODES 


IWEXTENSIONAL 
MODES 


0-1 0-3 0-4 0:5 0-6 


Fig. 8. 
Natural frequency in the extensional and inextensional modes. 
(Simply-supported edges, h/a=0-01) 
(1) Square plate, equal curvature in both directions. 
(2) Square plate, curved in one direction only. 
(3) Infinitely long rectangular plate, curved across the width. 


The shear strain arising from the bending is given by 

The lateral displacements u and v can be evaluated from the expressions for the 


direct strains and these displacements must be such that they are zero at the plate 
edges in the centre plane. They are 


wao- | ‘ (18) 


where f (9) is such that u=0 at ae and a/A 


end | ( 


and b/B 
z=0. 
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Here again, to obtain the frequency, the Rayleigh method is used and the mode 
considered for the case of rectangular plates, curved in both directions, with simply- 
supported edges, is 


w=J sin 


b 
sin (20) 


b 


which is asymmetric in both directions. 


This leads to lateral displacements such as 


The potential and kinetic energies can be evaluated as before and the natural 
frequency calculated, the main difference now being that the lateral displacements 
contribute to the total kinetic energy. For the particular case considered the 
natural frequency is given by 


a) 
(22) 


This expression reduces, for flat seal to the exact solution for the frequency 
of vibration when the mode is that of the first harmonic. Equation (22) is used to 
calculate the fundamental inextensional mode frequency only when the plate is 
curved in both directions. If the plate is curved in only one direction, then the 
fundamental inextensional mode is asymmetrical in the curved direction and 
symmetrical in the non-curved direction. All these expressions show that the 
frequencies fall slightly with increase in curvature. In Fig. 8, the natural frequency 
of the inextensional mode is compared with that of the extensional mode for the 
case of plates with simply-supported edges. It can be seen that for small curvatures 
the extensional mode is of lower frequency than the inextensional mode. As the 
curvature of the plates increases, the frequency of the extensional mode becomes 
higher than that of the inextensional mode. The inextensional mode falls slightly 
with curvature, the fall not becoming appreciable until a/A is about unity. 


Similar results would be obtained for plates with fixed edges, when a suitable 
expression for the inextensional mode, for plates curved in both directions, would be 


w = (b —2B¢) (a— AO)? (b — Bo)’. 
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A Note on the Evaluation of the Supersonic 
Downwash Integral 


B. A. HUNN, B.Sc. 
(lately of Hawker Aircraft Ltd.) 


Summary: To calculate tail loads in flight it is necessary to know the angle of down- 
wash at all points of the tailplane for any given value of wing incidence. Current 
trends in aircraft design place the tailplane off the plane of the wing. There exist 
solutions for the downwash in the z=0 plane for a delta wing with subsonic leading 
edges@: 4, This note gives a form of integral suitable for numerical evaluation which 
determines the downwash in the plane of symmetry (v=0) of a delta wing with subsonic 
leading edges. This note also points out an apparent error in Ref. 4 by G. N. Ward and 
gives a closed form for the downwash at a point whose forward Mach cone totally 
includes an arrowhead wing with supersonic trailing edges and whose centre line 
coincides with the axis of the cone. 

The procedure adopted is the use of standard elliptic function theory, which in 
certain cases can eliminate a large amount of tedious analysis. The problem has already 
been dealt with in general terms in Ref. 2, although the particular solution for the 

arrowhead wing is not included. 


Introduction 


In the determination of loads on an aircraft and the estimation of its response 
in manoeuvre it is necessary to know the aerodynamic loads on the wings and 
tailplane. The calculation of wing loads has been treated exhaustively by numerous 
authors. The determination of the aerodynamic influence of the wings on the 
tailplane loads, however, has not been dealt with to such an extent, largely because 
of the complexity of the integrals involved. 


The problem discussed here is that of downwash induced in the plane of 
symmetry of the aircraft when the wings have straight leading edges whose sweep- 
back is such that they are subsonic in a supersonic main stream. In the first part, 
the downwash due to a delta wing at incidence, having a simple triangular plan form, 
is considered. In the second part the case of swallow tail wings is discussed. 


Notation 
B? = M?-1 
E(k),E complete elliptic integrals of the second kind 
E’(\) = 
E(#,k),E(2,k) incomplete elliptic integrals of the second kind 
c maximum chord of the delta wing 
Cy,C, see after equation (29) 
d__ distance of point (x,0,z) behind trailing edge 
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complete elliptic integrals of the first kind 
incomplete elliptic integral of the first kind 
modulus of the elliptic functions 


| Jacobian elliptic functions 


Jacobi’s Theta function 

Jacobi’s Zeta function 

sn~' (iv/k) 

wing incidence 

see equation (37) 

see equation (31) 

tan y 

see after equation (29) 

semi-vertex angle of the wing 
= /z 

tan-' {(1/ Bz) [((x— x’)? — B*z?}!} 
= By 

velocity potential of disturbance 

z-component of the disturbance velocity 

Cartesian co-ordinates 

variables of integration 

K ({1 — k*}') 

see equation (19) 


finite part of the integral (see Ref. 1) 


2. Theory 


Using a right-handed system of Cartesian co-ordinates (x, y,z) such that x is 
measured positively in the direction of the undisturbed stream and z is positive 
downwards, with the origin at the vertex of the delta wing, the solution of the 


linearised lift problem which satisfies 


where Vo=u, the disturbance velocity vector, 


and B? = M? —-1, 


is given by 


on the wing, and the subsonic leading edge condition By < 1 is satisfied. 
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F (9, k) 
k 
sn (u, k) } 
cn (u, k) 
i dn (u, k) | 
(u) 
._Z(u) 
B 
y 
p 
w 
(x, y, 2) 
y’,u,x,t 
K’ (k) 
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SUPERSONIC DOWNWASH INTEGRAL 
It has been shown") that the velocity potential at any point inside the Mach 
cone from the origin is given by 


z(x— x) dx'dy’ 
[(y — y’}? +27] [(x — x’)? — B? (y — y’)? — B*z?}} 


=1((% 
— ay (3) 


where the integration extends only over the wing in the forward Mach cone from 
(x,y,z) where (2) holds, since in the vortex sheet of the wake 0¢/0x=0. To 
facilitate the analysis it will be assumed that ¢ (x,0,z) is required and furthermore 
that the delta wing is totally within the forward Mach cone from (x, 0, z), i.e. 


[(x —c)? — B*z?}! > Ac, 
where c is the maximum chord of the wing. The problem now is the evaluation of 


c nx’ 


(4) 


For this the following substitution is made: — 
= (u, k) 


snu being a Jacobian elliptic function. Writing »=»x’/z, we get 


c K 
2 z,U" | x’) dx’ k* du 


(x, 0, z)= 


mE’ (A) A (1+ 
0 


The integral in u can be evaluated’? and gives the form 


K K 

du | du 
| 
0 


where k*sn?a= —v?, 
This can be written”? finally as 


+ K(k) Z(2) } 


K (k)+ OKT 


cna 


where 9 and Z are Jacobi’s Theta and Zeta functions respectively. 
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However, dn (u, k)- ee 


and since dn uw and 9 (u) are even functions of u, (9) reduces to 


K (k) {14 


cna dna 
To evaluate the downwash velocity w, we use the result 


ao 
0z 


that | 


K (k) k (x— x’) 
zA 


*_7(2)} dx. 


Splitting this integral into two parts and using the definition of Z(z) and sn z gives 


c 


2 
2 2,U1? Kw kav : (14) 


22 
0 


c 
2 ix’ 


(KE(2,k)—2E} dx’. (15) 


using the definition of KZ (z) and the fact that k sn z= iv. 


To perform the differentiation, the following results are easily obtainable : — 


i(x—x’)k 


_ k’sna =) 
dnz 


k)—(1—k*)z 


K} 
ex? (1—k*) 


zk? (K —E) 
nex’? 
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Substituting these results where necessary in (14) and (15) and adding, the final 
result is obtained :— 


2a,U x’) 1 , 


(KEG, k) dx’ } 


However, 


i (KE(, k) aE) —KF KAx’ 


where o=tan—’ 


[(x — x’)? — 


Thus on combining equations (16), (17), (18), and using standard tables of 
incomplete elliptic integrals it is now possible to evaluate the integrands at as many 
values of x’ as are necessary to obtain a good approximation to the integral. 


It is interesting at this stage to pursue the limiting process of letting z tend to 
zero. Writing 


we consider 
lim | { x2 4 KAx”? 


E(x-—x)k( k? } 


lim k (x — x’/)=Ax’ 


lim = = KE’ + EK’ — KK’=2/2. 


z 0 


c 


Ez z d= 
0 0 


and, since & is a finite continuous function of x’ within the range, it may be shown 
that the limit of the integral is simply 


: [n?x’? + 
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11) 
12) 
13) 
Now 
} But 
|_| 


w (x, 0,0)= { 


Ax’ xr 


But since = Gx)’ 


22U {(A(K-E) 


Ww 0, 0) mE’ (A) 


Ac/d 

_ 24U f 

= + dk}, (27) 
0 


These are the results obtained in Ref. 3 and quoted finally in Ref. 4. The results 
obtained in Ref. 4 equation (12), will be seen to be incorrect, which is apparently 
due to a fallacious statement in equation (11). 


Consider the downwash due to a swept wing, when the pressure on the lower 
surface is given, using 


2,Un? x 
Ox TE’ [n?x? — y?}}! 


(28) 


In order that this may be true we have to make the following restrictions on 
the wing: 


(i) The trailing edges must be supersonic with respect to the trailing edge 
apex Mach cone for the speed range considered. 


(ii) The tip must be pointed, so that in fact swallow tail wings are considered. 


Now the downwash will be obtained from the integral as shown in (3), with 
the bounds of integration modified. However, it is equally simple to complete the 
delta by filling in the triangular wake region whose vertices are the two wing tips 
and the trailing edge apex and, using the downwash formula for complete delta 
wings, to obtain that due to the swallow tail wing by subtracting the downwash due 
to this triangular wake as if it were wing. 


Thus if we require the downwash at A (x, 0, z) behind the swallow tail we have 
to evaluate the correction term 


n,(@’ 


(29) 
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where x=c, defines the trailing edge apex, 
x=c, defines the value of x at the tips, 
and =cot Ay, 


the sweepback of the trailing edge being A;. There is one other condition which 
has to be accounted for, namely that involving the amount of the wing included in 
the forward Mach cone from A. We shall consider here the algebraically more 
simple case when the whole wing can be included. This can be expressed by 
the condition 


yx’ <(1/B)(x—x’), when z=0. . (30) 

We use the same substitution as before and introduce £ by the definition 
yx’ sn B=, (x’—c,). G1) 
The general solution in the xz-plane has been obtained, but owing to. its algebraic 


complexity it has not been included. The correction term for the downwash at 
A (x,0,0) due to an arrowhead wing is found to be, using a result from Ref. 1, 


U rE’(A) x’sn* (u, x—x’ 


w’ 2 2, | # dx'du 
Ss 


where S’ denotes the triangular cut-out, and .Y indicates “finite part.” We have :— 


B 


du 

| F (¢, k)—-E(¢, k) 
0 


cnB dng 


(33) 


/ 


x x’snB 


x-Cc 
where snp = 
nx 


Hence the final result for the arrowhead wing totally included in the forward Mach 
cone from A is 


w 2a {[(K-E),, 


& 


+ 


c 


xn, (x’—¢,) 
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The last integral whose principal part is required can be evaluated by the standard 
process of writing 


where 2’, y’, 5 are the roots of oF —X,? (x - in x. 
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The Design of Thin Finite Wings in 
Incompressible Flow 


G. J. HANCOCK, B.Sc. 
(Department of Mathematics, Sheffield University) 


Summary: A method is developed for designing thin finite wings with prescribed 
lift distributions. Particular attention is given to wings which, at a certain optimum 
angle of incidence, have zero loading along the whole leading edge. It is suggested 
that such a wing will have a greater range of lift coefficient for low drag than other 
wings with the same plan form. 

On the basis of the linearised theory a formula, relating the wing surface to the 
pressure distribution over the surface, is deduced which differs from the usual potential 
relationship by avoiding the difficulty of the non-uniform convergence with respect to 
the normal distance away from the wing surface singularities. The surface integral 
which is developed includes a regular singularity and a definite limiting process. 

Three examples have been worked out in which the condition of smooth flow at 
the leading edge is satisfied. 

In the first example, the pressure distribution of a rectangular aerofoil is assumed to 
be elliptic in the chordwise direction and parabolic in the spanwise direction. The 
reverse pressure distribution (parabolic chordwise, elliptic spanwise) is considered in the 
second example over the same rectangular aerofoil. Thirdly a pressure distribution 

is prescribed over a swept-back wing. 


Introduction 

Of the two basic aims of aerodynamic calculation, namely 

(i) to determine the pressure distribution in flow past an obstacle of given 

shape 

(ii) to determine an obstacle shape having a given pressure distribution, 
only problem (i) has received great attention in problems of flow over three- 
dimensional wings at low Mach numbers, although both problems have been 
extensively treated for two-dimensional aerofoils. 


As is well known, the symmetrical part of the pressure distribution, resulting 
from the fairing of the wing, is much easier to determine than the anti-symmetrical 
part, resulting from camber and angle of incidence. This is because, to the first- 
order approximation used, the symmetrical part of the potential functions is deter- 
mined from a boundary condition which is of the same type all over the plane which 
contains the wing; while the boundary condition for the anti-symmetrical part is of 
mixed type. This is unfortunate because the lift distribution is what is most needed. 


For problem (ii) the situation is reversed. The easier problem here is the 
determination of the shape of the central surface such that the lift at a given angle 
of incidence is distributed in a given manner. (This is the extension to three 
dimensions of Goldstein’s theory of the design of camber lines.) 
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No work seems to have been done on problem (ii), while theoretical work on 
lift distribution on wings has consisted either of accurate computations taking many 
months, or approximate methods with restricted ranges of application. The labour 
of the accurate methods has been justified as helping to indicate more precisely the 
ranges of application of the approximate methods. 


Now if, instead, some exact solutions are found by solving the design problem 
(ii), with (as will be seen) far less labour, these can equally be used as a testing 
ground for approximate methods of the solution of problem (i). This in itself 
makes the development of the necessary technique worth while. In addition, 
theoretical knowledge of what sort of camber distribution and twist may be expected 
to produce specific effects on lift distribution, particularly for unconventional plan 
forms, is almost certainly worth having. 


The question now arises—what sort of lift distribution should be aimed at? 
In this first paper on the subject the author has decided to concentrate on designing 
a central surface such that each section of the wing achieves its optimum incidence 
simultaneously, namely at the cruising lift coefficient. If a suitable fairing were 
added, this wing would have quite a range of values of C, about this optimum C,, 
in which the drag might be expected to be low. By optimum incidence is meant the 
condition of zero loading at the leading edge so that there would be no tendency to 
breakdown even if no fairing whatever were put onto the central surface. 


The general problem of determining the shape of the central surface from the 
lift distribution is solved in Sections 2 and 3. In this work the fairing is neglected 
and the problem reduces to the determination of the shape of an (infinitely) thin 
wing from its assumed pressure distribution. The difficulty of the simplest formal 
approach, in leading to a result containing improper integrals, is avoided by an 
extension of the method devised by Ursell® for the determination of the pressure 
distribution on symmetrical wings at zero lift. As a result the determination of a 
wing surface for a prescribed pressure distribution becomes a direct double 
integration. 


In Sections 4, 5 and 6 three examples satisfying the condition mentioned, of 
zero wing loading at the leading edge for a prescribed value of C,, are designed. 


Notation 
x,y,z rectangular Cartesian co-ordinates 
z=((x,y) equation of the wing surface 
x=I(y) | equations of the leading and trailing edges of the wing surface on 
x=t(y) | the plane z=0 
V __ velocity of the main stream 
(u,v,w) disturbance velocity at point (x, y, Z) 
p(x,y,z) disturbance pressure 
(x,y,z) velocity potential 
Sw wing projection on plane z=0 
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S, trailing vortex projection on plane z=0 
R_ vector X, Y,0 
r vector x, y, Z 
Sp small surface in space z > 0 around (X, Y, 0) 
Cp curve of intersection of S, with xy-plane 
S;, small surface in space z > 0 around x=/(y) 
S,, small surface in space z > 0 around y= +s 
S_, small surface in space z > 0 around y= —s 
S’ large hemisphere in space z>0 centre (X, Y,0) 


Sw’+Sy’ — the plane surface Sy +S; excluding the interior of curve C; and the 
regions of the perimeter, x=/(y) and | y|=s 


¢,a__ sides of the rectangle C; parallel to x- and y-axes respectively 
P, limiting operation 0 first and a—>0 secondly 


a angle of incidence 


2. Mathematical Formulation 


The incompressible flow past a thin wing at a small angle of incidence is now 
considered. As shown in Fig. 1, the origin of the rectangular Cartesian axes is 
taken at the mid-point of the centre chord of the wing. The equation of the wing 
surface is defined as z=((x, y), and the equations of the leading and trailing edges 
of the projection of the wing surface on the plane z=0 are x=/(y) and x=1?(y) 
respectively. 


If V is the velocity of the main stream at infinity, and (u,v, w) is the disturbance 


(y) 


n=tly) 
| 


Fig. 1. 
Plan and elevation of a thin finite wing, 
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velocity at the point (x,y,z) due to the presence of the aerofoil in the stream, the 
linearised equations of motion are 


ou __ 

pox 

op 
p oy (1) 

Ox p 0z 


The equation of continuity becomes 


Ou dw 


The problem is to determine the relationship between a doublet distribution over 
the projection of the wing and trailing vortex sheet on the plane z=0, and the shape 
of the wing surface. The conditions to be satisfied are that smooth flow exists over 
the wing surface and that the pressure is continuous over the vortex sheet. The 
condition of adjacent flow states that on the wing plan form 


w (x, y, 0) -— a (x, y) 


where w (x, y,0) is calculated from the doublet distribution. 


The flow outside the region of the doublet distribution is irrotational, therefore 
a velocity potential ¢ exists such that 


Then the problem reduces to the determination of the discontinuity in ¢ over a 
certain surface in the xy-plane with the condition that (0¢/0z),-, satisfies 
equation (3). 


Denoting the wing projection on the xy-plane by Sw and the projection on the 
xy-plane of the trailing strip from x=1(y) to x= + 0, (| y| <s) by Sy, the following 
boundary conditions are to be satisfied. 


(i) Since ¢ is an odd function of z, i.e. o(x, y,z)= —9(x, y, —2), 
¢=u=v=0. 
w remains finite on the plane z=0 outside Sy + Sr. 


(ii) Over Sw; the discontinuity ¢,-,,—¢%.--, is a function of x and y; therefore 
discontinuities in u,v exist. Since 0¢/0z is an even function of x, w is 
finite, satisfying equation (3). 
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(iii) Over Sy, since there is no pressure discontinuity across S;, the velocity 
u must be continuous and therefore u=0, hence ¢,_,, —¢%.-_, is a function 
of y only. 


(iv) It is assumed in accordance with the usual linearised theory that in the 
region of the perimeters of Sw + Sr, 


— = O (54) 


when 6 is the distance from the particular perimeter, either near the 
leading edge x=/(y) or the lines | y|=s. 


Thus Uza49 — =O (6-4), Vz=49 — Vz= ~=O6 ~4) 


a 
and “Uz=-y) =O (8 dy ) 


near the boundaries mentioned. 


3. Application of Green’s Theorem 


Following Ursell’ a point P(X, Y) on Sw is considered. The point (X, Y,0) 
is denoted by the vector R and (x,y,z) by r. A small surface S$; is constructed 
around P in the space z>0, cutting the plane z= +0 in a curve Cp. Green’s 
theorem is to be applied to the two harmonic functions w and 1/(R—r) over the 
region bounded by the plane z= +0 (excluding singularities) and a large hemisphere 
above z=0. 


Since —— = — — + — 3/2), 


from boundary condition (4), in the region of any perimeters of Sw+Sy, it is 
necessary to construct small cylindrical surfaces S,, S,, S_, around these perimeters 
as shown in Fig. 2. 


If S’ is a large hemisphere, centre (X, Y, 0), 


where the integration is over the surface comprising S’, Sp, S,, Sy., S_,, and the 
plane z= +0 exterior to C, which completes the closed surface. If the radius of 
S’ tends to infinity, the integral over S’ vanishes and equation (4) reduces to 


the integral being taken over the surface S,+5,,+S_,+Sp+Sw’'+Sy where 
Sw’ +S,’ denotes the surface Sy +S, for z= +0, excluding the interior of Cp and 
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Fig. 2. 


Surfaces around the singularities. 


excluding the regions of the perimeters. Each of the integrals of (5) will now be 
considered separately. 


(a) First the integral over z= +0 will be considered. 


since 0/d0n= —0/0z. 
Replacing 0w/dz by —(0u/0x+0v/d0y) and integrating by parts, 


{| 1 ow 
) 
+S,’ 


+8 


[(x— X)? +(y— Y)*]}! (x XP YFP" 


P we 


dxdy . (6) 


c 


where dx/d£=sin 6, dy/d&=cos@ around the curve 


(b) Secondly the integral over S; is considered, where Cp is taken as a rectangle 
of sides 2<, 2a parallel to the x- and y-axes respectively. S, is now taken as a 
semi-cylinder of length 2a, radius «. The integral over the surface S, then reduces 
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to 2xw (X, Y) if the limiting process : —> 0 is applied first and the process a—> 0 
afterwards. 


(c) Finally the integrals over S,, S_,, S,, are evaluated, as these surfaces tend 
to zero and may be shown to cancel out the line integrals over the perimeter of 
Sw+ Sp in equation (6). 


Hence from (a), (b) and (c), equation (5) reduces to 


(x, y)cosé+v(x,y)siné ,. 


j 


OY } 


SwtSq 


where P;, denotes the operation lim « —> 0 first and lima—> 0 afterwards. 


Since the line integral vanishes in the application of P,, 


(x—X) u(x, y)+—Y) v(x, y) 
22w(X, Y)- { | dxdy \ 
SwtSp 


4. Determination of a Wing Surface of a Rectangular Wing of 
Aspect Ratio 2 with Optimum Angle of Incidence along 
the Span 


It was explained in Section 1 that it is desirable to design a wing possessing 
smooth flow along the leading edge. This is performed by prescribing a pressure 
distribution which vanishes like 5! (where 4 is the small distance from the edge) at 
the leading edge as well as the trailing edge. By the relationship of equation (1), 
the velocities u and v on the surface Sw + S; are determined; therefore on substitution 
in the right hand side of equation (7) the velocity w (X, Y) is determined, and hence 
¢ (x, y) is known by equation (3). 


To simplify some initial calculations, the case of a rectangular wing is con- 
sidered, where the plan form is given by -l1<=x< +1, -s<y<+s. 


The pressure distribution is assumed to be 
V*Kz(1- x)! (1 
a (1 — x’) 


where K is a constant such that the angle of incidence at the centre section is z. 
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Hence by equation (1), 
u(x,y) _ y? 
=Ka(1-x(1-¥) on Sw 


=() on Sy 


_ Al | 
x (1 —x’}!+sin on Sw 


=-Ke-= on Sy. 
Ss 


On substitution in equation (7), 


(X,Y) 

P { X) (1 — (1 — y?/s?)—(y/s*)(y — Y) [x (1 — + sin-' 
L [(x- XY +(y- 


Sp 


In this formula s=2 is taken. The integration with respect to y is performed; then 
the integration with respect to x is calculated numerically. The values of 
w(X, Y)/(VKz2) are given in Table I, and the graph of the surface z=((x, y) is 
shown in Fig. 3. It must be remembered that equation (3) gives only the relationship 
between 0¢/0x and w(x, y,0) due to the assumptions of the linearised theory, and 
therefore in the determination of ((x,y) an arbitrary function of y has to be 
assumed. In this example the assumption is made that the leading edge lies on the 
xy-plane. 


From the graph in Fig. 3 the value of K is 3-6. 


TABLE I 
VALUES OF w(X, Y)/(VKa2) 
00 0:5 1:0 2:0 
X= -—-1:00 0-673 0-637 0-625 0-407 0-147 
—0-75 0-436 0-434 0-359 0-318 0-083 
—0-50 0-178 0-177 0-160 0-201 0-285 
—0-25 -0:063 — 0-054 — 0-006 0-107 0-324 
—0-00 — 0-289 — 0-262 -0-164 0-042 0-348 
+0-25 — 0-468 — 0-434 — 0-284 0-007 0°426 
+0-50 — 0-643 — 0-590 — 0-386 — 0-009 0-521 
+0-75 — 0-788 —0-714 —0-476 —0-021 0-625 
— 0-760 —0-701 
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Fig. 3. 
Wing Surface z=( (x, y). 
Rectangular wing with optimum incidence across the span. 
n 
yf 
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Fig. 4. 


Wing surface z=( (x, y). 
Rectangular wing with elliptic circulation distribution. 
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p= —pV?K2(1 ): 


Hence “(1 on 
=0 on 

V~ #0-y 373 


(1 — y?/s?)! 


on 


K=2°7. 


be applied. 


TABLE Il 


VALUES OF w/(V Kz) 


5. Determination of a Rectangular Wing of Aspect Ratio 2 with 
Prescribed Elliptic Circulation Distribution 


In the last example a pressure distribution is prescribed which satisfies the 
flow conditions in the chordwise direction but which states that the circulation 
distribution in the spanwise direction is parabolic. The wing which satisfies this 
pressure distribution is seen to have a total twist in the region of the wing tips. To 
show that this twist is due to the assumption of the parabolic distribution of 
circulation, an example is now worked out for an elliptic circulation distribution, 
together with a parabolic chordwise distribution to simplify the calculations. 


The prescribed pressure distribution is therefore taken to be 


Sw 


Sy 


Sw 


Sy. 


The values of w/(VKz) have been computed at the two sections Y=0, Y=1-75 
and are given in Table II. The graph of this wing is indicated in Fig. 4, where 


It is seen that, to achieve a uniform pressure distribution in the spanwise 
direction, no twist is needed and only a slight reduction of camber at the tips should 


A=f0 - 1-60 — 1-40 
0-5 —0-80 —0-70 
0-0 —0-33 —0-33 
—0°5 —0-05 +0-05 
—1:0 +0-90 +0-70 
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Delta wing. 
6. Determination of a 45° Swept-back Wing 
The greatest interest is in the application of the design criterion to swept-back 
wings. One example is now worked out, the delta wing whose equations of the 
leading and trailing edges are 
x= -6+|y| 
and x-1l= 0, 
15 and whose span s=6. This wing is shown in Fig. 5. 
ere The prescribed pressure distribution will be assumed to possess the properties 
previously stated (i.e. elliptic circulation, smooth flow at the leading and trailing 
f edges), together with the property that the derivatives of the pressure across the 
id centre line of the aerofoil are continuous. Therefore the pressure distribution is 
. taken to be 
p=pV* x +64 5 1 ) 
TABLE III 
VALUES OF w/(VKz) 
Y= 0 2 4 
X=-6 + 0-04 
—0-03 +0-10 
X=-2 —0-18 —0:01 +0-15 
X=-1 -0-50 —0-25 +0-05 
—0-80 
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Ys 


-2 


Fig. 6. 
Wing surface z=(¢(x, y). Delta wing. 


Therefore 


on Sw 


on S; 
Vv yKz2 1 2 | 27 


x [1 —x) (x+6—|y|)#- } on Sw 


— 


The values of w/(VKa) have been computed at the three sections Y=0, 2, 4 and 
are given in Table III. In the graph of the wing (Fig. 6) the assumption has been 
made that the trailing edge lies on the xy-plane, and it is seen that K=2°8. 

The wing seems to have only a slight constant twist in the spanwise direction, 
as the camber at each section increases. It will be noticed that the leading point 
of the delta wing is hardly curved at all to meet the oncoming fluid. 
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